Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  796 1728 (N59E03 L=  2)  1988    301        4 ( 706/1650- 914/1788)  10.9 E14-W12 N53-N62
      2 1242 1721 (N46W38 L= 43)  1131    183        0 (1152/1690-1329/1748)  -9.5 W29-W47 N43-N48
      3 1466 1692 (N38W56 L= 60)   746     92        2 (1426/1666-1511/1711) -20.6 W50-W61 N35-N39
      4 1561 1520 (N25W50 L= 55)   365    131        1 (1523/1483-1596/1545) -36.9 W45-W52 N23-N27
      5  649 1345 (N32E20 L=345)   258     77        1 ( 620/1314- 688/1364)  11.3 E22-E17 N30-N33
      6 1379 1262 (N14W27 L= 31)  1085    252        5 (1327/1158-1414/1355)  80.5 W23-W29 N08-N20
      7  632 1090 (N17E24 L=341)   420    151        0 ( 563/1063- 698/1120)   5.7 E28-E18 N15-N17
      8  508  648 (S08E40 L=324)   305     72        0 ( 468/ 628- 550/ 665) -18.4 E43-E37 S09-S06
      9  966  225 (S48E26 L=339)   488     29        0 ( 944/ 206- 991/ 249) -90.0 E28-E23 S49-S44
  

    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