Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  485 1577 (N15E49 L=227)   357     38        0 ( 461/1551- 510/1597)  56.3 E50-E46 N13-N16
      2 1219 1414 (N20W03 L=279)   991    221        5 (1135/1333-1299/1485) -20.0 E03-W08 N15-N23
      3 1227 1183 (N06W08 L=284)   284     82        1 (1188/1146-1261/1220) -31.0 W05-W10 N04-N07
      4 1761  979 (N08W48 L=324)   653    174        1 (1728/ 922-1817/1012)  71.6 W45-W53 N04-N10
      5  248  895 (S29E52 L=224)   343     96        1 ( 221/ 871- 281/ 921)  21.8 E54-E48 S30-S26
      6  314  817 (S33E45 L=230)   392    106        0 ( 272/ 804- 355/ 829)  14.0 E49-E41 S33-S31
      7 1728  539 (S16W62 L=338)   399     64        0 (1697/ 508-1749/ 568)  78.7 W60-W63 S18-S13
      8 1474  467 (S27W45 L=321)  1135    262        0 (1352/ 436-1592/ 505)  35.5 W38-W52 S32-S22
      9  860  209 (S66W19 L=295)   514     30        0 ( 839/ 183- 881/ 231) -51.3 W15-W23 S68-S63
  

    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