Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1401 1670 (N31W47 L=262)   776    185        1 (1307/1648-1487/1691) -24.8 W39-W52 N27-N34
      2  879 1493 (N35W03 L=218)  1771    633        0 ( 607/1432-1161/1541) -14.7 E17-W21 N27-N38
      3 1445 1337 (N11W33 L=248)   848    208        0 (1391/1243-1490/1420) -90.0 W31-W34 N04-N16
      4 1863 1168 (S10W62 L=277)   570     78        0 (1844/1133-1883/1211) -90.0 W60-W64 S11-S07
      5  355  879 (N10E46 L=169)  1761    513        0 ( 209/ 695- 523/ 996) -46.4 E58-E37 S03-N18
      6 1718  623 (S39W42 L=257)  2899    454        4 (1584/ 493-1831/ 812)  30.0 W27-W53 S45-S29
      7 1350  713 (S22W12 L=227)   780    220        1 (1315/ 628-1381/ 800)  52.1 W08-W15 S26-S17
      8  506  403 (S21E54 L=161)  1126    195        5 ( 477/ 320- 547/ 496)  47.5 E59-E46 S27-S15
      9  659  420 (S23E41 L=174)   409    151        0 ( 641/ 343- 684/ 478)  41.6 E45-E35 S28-S20
     10 1218  278 (S50E10 L=205)  1120    205        4 (1131/ 236-1309/ 326)   0.0 E18-W00 S52-S47
  

    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