Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1112 1863 (N55E09 L=342)   532     84        0 (1085/1836-1162/1882)  39.8 E11-E04 N51-N58
      2  458 1707 (N32E59 L=292)   495     63        0 ( 445/1669- 473/1741) -36.9 E63-E55 N29-N35
      3 1188 1505 (N26W04 L=355)   488    154        0 (1156/1441-1211/1575) -71.6 W01-W05 N21-N30
      4 1792 1252 (N20W51 L= 42)  1005    230        1 (1756/1143-1820/1353) -71.6 W48-W54 N13-N25
      5  404  951 (S16E40 L=310)   269     56        1 ( 387/ 917- 426/ 979)  53.1 E41-E38 S18-S14
      6  456  845 (S22E36 L=315)   298     95        0 ( 421/ 813- 472/ 889)  59.0 E38-E34 S24-S19
      7  627  791 (S24E23 L=328)   808    245        1 ( 558/ 676- 662/ 889)  66.8 E28-E20 S32-S18
      8 1118  195 (S61W34 L= 24)   567     57        2 (1097/ 166-1139/ 225) -56.3 W29-W39 S63-S57
  

    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