Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  739 1741 (N62W05 L=244)   568    101        1 ( 704/1702- 781/1778)  14.0 W00-W12 N58-N64
      2 1392 1700 (N34W52 L=290)  1667    290        4 (1264/1652-1528/1744) -38.7 W43-W60 N26-N40
      3 1109 1642 (N40W28 L=266)   332     64        0 (1089/1611-1139/1675)  78.7 W26-W29 N37-N42
      4  455 1199 (N31E33 L=205)   395    113        1 ( 399/1177- 507/1218)  -5.7 E38-E28 N29-N31
      5 1840  918 (S23W54 L=293)   642     77        1 (1821/ 876-1861/ 953)  45.0 W51-W56 S25-S21
      6 1797  801 (S29W48 L=287)   997    143        1 (1760/ 753-1862/ 841)  -4.4 W43-W57 S32-S25
      7 1129  657 (S17E04 L=235)   383    109        0 (1083/ 639-1187/ 676) -18.4 E06-W00 S17-S14
      8 1075  504 (S25E11 L=227)   670    270        1 (1008/ 456-1131/ 567) -29.7 E15-E08 S27-S20
      9  649  386 (S22E44 L=195)  1251    190        3 ( 578/ 355- 733/ 435) -14.0 E49-E37 S25-S17
  

    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