Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  382 1620 (N20E64 L=336)   647     51        1 ( 359/1600- 408/1641)  11.3 E66-E60 N18-N21
      2  459 1493 (N13E49 L=350)  1509    397        3 ( 402/1388- 544/1607) -32.3 E59-E38 N08-N20
      3 1007 1408 (N20E10 L= 30)   611    169        0 ( 962/1346-1059/1475) -29.1 E14-E05 N16-N23
      4 1430 1339 (N26W18 L= 58)   329    106        1 (1397/1284-1460/1381) -39.8 W15-W21 N22-N27
      5 1032 1216 (N09E04 L= 36)   366     76        3 (1004/1179-1066/1253) -36.9 E05-W00 N07-N10
      6  413 1203 (S05E43 L=357)  1009    181        3 ( 337/1151- 491/1244)  36.0 E48-E37 S08-N00
      7 1748 1163 (N23W46 L= 85)   540    138        0 (1709/1120-1782/1188) -18.4 W42-W49 N20-N24
      8  658 1030 (S10E22 L= 18)   637    190        1 ( 576/ 989- 749/1066)  10.3 E27-E16 S11-S07
      9 1774 1007 (N14W50 L= 90)   491    116        3 (1744/ 949-1794/1045) -71.6 W47-W53 N10-N16
     10 1361  918 (S01W22 L= 62)   254    106        1 (1311/ 882-1407/ 946)   0.0 W18-W25 S02-N00
     11 1130  602 (S25W17 L= 56)   350    116        0 (1073/ 585-1189/ 616)  15.9 W13-W20 S26-S23
     12  676  555 (S39E13 L= 27)   235     42        0 ( 651/ 540- 707/ 570)  14.0 E14-E10 S40-S37
  

    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