Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1339 1374 (N24W10 L= 36)   726    237        3 (1282/1277-1399/1479) -36.9 W03-W15 N19-N29
      2  887 1282 (N07E14 L= 12)   329    130        0 ( 834/1255- 965/1313)  45.0 E16-E10 N04-N10
      3 1188  980 (S02W10 L= 37)   185    113        0 (1141/ 962-1242/1022)   0.0 W06-W13 S03-S01
      4  528  884 (S24E27 L=360)   763    210        0 ( 448/ 847- 607/ 951)  45.0 E33-E22 S27-S18
      5 1159  816 (S12W13 L= 40)   256    124        1 (1111/ 795-1229/ 842)  45.0 W10-W16 S14-S08
      6 1634  684 (S06W48 L= 75)   376     60        1 (1607/ 653-1651/ 713)  90.0 W45-W49 S08-S04
      7  581  617 (S40E19 L=  8)   996    179        2 ( 511/ 537- 656/ 682) -23.2 E25-E11 S43-S35
      8 1481  588 (S16W41 L= 68)   400    134        0 (1457/ 531-1527/ 659) -71.6 W39-W43 S19-S10
      9  823  549 (S37W01 L= 28)   323     66        1 ( 781/ 528- 864/ 567)   9.5 E02-W04 S38-S36
     10 1321  530 (S24W33 L= 60)   707    275        1 (1207/ 512-1458/ 547)  23.6 W25-W41 S27-S19
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y