Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  867 1545 (N32E18 L=287)   388     93        1 ( 840/1506- 911/1574)  53.1 E20-E14 N29-N34
      2 1803 1378 (N32W63 L=  8)   496     54        1 (1783/1352-1823/1406) -56.3 W59-W65 N30-N33
      3 1072 1234 (N14W01 L=306)   404    184        2 (1018/1165-1101/1324) -60.9 E03-W03 N10-N19
      4 1487 1114 (N11W29 L=334)   409    209        0 (1454/1018-1511/1208) -70.0 W25-W31 N05-N16
      5 1735 1080 (N12W49 L=354)  1209    247        4 (1715/ 965-1754/1200)  82.4 W47-W51 N04-N19
      6  589  935 (S09E27 L=279)   462    124        1 ( 565/ 880- 613/ 995)  81.9 E27-E25 S12-S05
      7  188  852 (S19E64 L=241)  2783    227       12 ( 142/ 768- 233/ 943) -90.0 E75-E57 S24-S13
      8  307  847 (S18E50 L=256)   356     68        2 ( 271/ 829- 333/ 867)   0.0 E53-E46 S19-S16
      9 1640  692 (S13W47 L=352)   312     53        1 (1616/ 664-1658/ 720)  56.3 W45-W47 S14-S11
     10 1381  598 (S22W29 L=335)   271    102        1 (1344/ 571-1412/ 646)  68.2 W27-W31 S23-S18
     11 1692  476 (S27W64 L= 10)  3846    337        6 (1560/ 400-1805/ 576)  37.9 W52-W76 S33-S18
     12 1076  456 (S35W11 L=317)   273     88        1 (1063/ 418-1095/ 497) -59.0 W10-W13 S37-S32
  

    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