Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  787 1706 (N57W05 L= 42)   405     52        2 ( 763/1679- 812/1728)  26.6 W02-W08 N55-N59
      2 1202 1356 (N21W19 L= 56)   198     67        0 (1162/1339-1235/1369) -36.9 W17-W21 N19-N23
      3  958 1162 (N17E01 L= 36)   723    201        1 ( 882/1106-1020/1236) -68.2 E03-W02 N12-N22
      4 1784 1166 (S06W54 L= 91)   339     45        1 (1771/1140-1804/1194)  45.0 W52-W55 S07-S03
      5 1018 1001 (N06E01 L= 36)   464    130        1 ( 986/ 940-1043/1064)  56.3 E03-W01 N03-N09
      6  669 1009 (N14E21 L= 16)   213     98        0 ( 656/ 971- 694/1055)  53.1 E22-E18 N12-N16
      7  299  978 (N19E51 L=346)   359     70        1 ( 285/ 932- 309/1014) -78.7 E51-E49 N16-N21
      8 1765  904 (S20W46 L= 83)  1402    251        1 (1711/ 778-1818/1002)  42.5 W39-W53 S26-S15
      9  234  974 (N20E58 L=339)   347     45        0 ( 217/ 949- 253/1000)  45.0 E60-E55 N18-N21
     10 1867  896 (S25W59 L= 95)   848    119        0 (1826/ 846-1908/ 941)   8.7 W52-W65 S26-S22
     11  920  827 (S02E11 L= 26)   389    100        2 ( 898/ 776- 939/ 880)  68.2 E11-E08 S04-N01
     12 1770  758 (S30W46 L= 83)   489     91        0 (1731/ 741-1803/ 782) -14.0 W41-W49 S31-S27
     13 1128  611 (S19E04 L= 33)   419    168        0 (1054/ 582-1191/ 630) -26.6 E09-W00 S22-S17
     14 1326  412 (S39W03 L= 40)  1353    255        5 (1225/ 329-1397/ 496)  26.6 E06-W10 S43-S34
     15 1058  432 (S29E14 L= 23)   292     89        1 (1020/ 413-1104/ 449)  -9.5 E16-E10 S30-S27
     16  411  388 (S19E66 L=331)   591     73        2 ( 386/ 366- 432/ 410) -76.0 E70-E62 S21-S17
     17 1142  294 (S43E14 L= 22)   367     48        1 (1121/ 268-1165/ 324) -90.0 E16-E13 S45-S39
  

    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