Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1146 1726 (N52W11 L=211)   429     72        1 (1113/1696-1188/1746)   7.1 W07-W15 N49-N54
      2  513 1530 (N35E43 L=157)  1331    261        2 ( 413/1489- 628/1591)  -6.6 E57-E31 N32-N38
      3  722 1473 (N32E23 L=177)   328     97        0 ( 682/1445- 754/1512) -35.5 E27-E20 N29-N34
      4  981 1142 (N11E03 L=197)   186    102        0 ( 961/1098-1007/1181) -71.6 E04-E01 N07-N13
      5  778  685 (S18E16 L=184)   499    215        1 ( 680/ 663- 873/ 722)  -8.7 E22-E09 S19-S16
      6  417  569 (S28E47 L=153)  2258    337        4 ( 332/ 413- 498/ 686)  28.2 E77-E36 S41-S19
      7  833  610 (S23E12 L=187)   359    174        0 ( 781/ 548- 904/ 651)  37.9 E16-E07 S27-S20
      8 1593  423 (S37W51 L=251)   662     82        2 (1570/ 384-1613/ 460) -36.9 W46-W56 S40-S34
      9 1476  362 (S42W42 L=242)  1205    131        2 (1414/ 333-1541/ 390)  13.0 W36-W49 S44-S39
     10 1251  310 (S46W22 L=222)  1207    215        2 (1161/ 288-1335/ 344)   3.4 W13-W30 S48-S43
  

    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