Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  629 1669 (N60E12 L=280)   404     56        0 ( 613/1636- 643/1702)  81.9 E14-E09 N56-N63
      2 1599 1613 (N22W62 L=354)   524     87        0 (1562/1593-1642/1639) -51.3 W59-W64 N19-N24
      3  334 1561 (N59E59 L=233)   846     79        2 ( 286/1540- 371/1583)   0.0 E75-E49 N57-N60
      4  547 1404 (N42E24 L=269)   588    146        1 ( 519/1350- 614/1455)  24.0 E26-E16 N38-N43
      5  857 1357 (N30E01 L=291)   680    130        2 ( 808/1316- 918/1405) -66.8 E04-W02 N27-N34
      6 1844  960 (S21W55 L=347)   505     64        1 (1818/ 934-1876/ 985)   0.0 W51-W59 S22-S19
      7 1206  472 (S31E04 L=288)   358    123        1 (1178/ 439-1274/ 505)   8.1 E06-W01 S32-S29
  

    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