Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  657 1575 (N35E24 L=325)   300     45        1 ( 632/1554- 675/1601)  45.0 E25-E22 N33-N36
      2  633 1412 (N24E24 L=325)   228     70        0 ( 617/1380- 654/1450) -68.2 E25-E22 N21-N26
      3 1671 1360 (N15W46 L= 35)   399    102        1 (1636/1330-1709/1402) -51.3 W43-W48 N13-N18
      4 1459 1267 (N10W29 L= 17)   335     59        1 (1439/1236-1487/1303) -76.0 W27-W30 N08-N12
      5 1020 1174 (N07W00 L=349)   935    247        2 ( 912/1078-1101/1252) -50.7 E05-W04 N00-N11
      6  614 1120 (N06E25 L=324)   207     58        0 ( 582/1099- 647/1140) -26.6 E26-E22 N04-N06
      7  633 1057 (N02E24 L=325)   337    157        2 ( 602/1036- 675/1092) -26.6 E25-E21 N00-N03
      8  187  866 (S06E62 L=287)   660    148        2 ( 154/ 825- 224/ 900)  16.7 E67-E57 S08-S04
      9 1638  820 (S17W40 L= 29)   347    150        2 (1611/ 780-1655/ 869)  68.2 W37-W41 S19-S14
     10 1687  577 (S33W48 L= 37)   382     72        1 (1663/ 550-1720/ 593)  36.9 W45-W51 S34-S31
     11 1410  442 (S41W27 L= 16)  1618    424        7 (1224/ 372-1592/ 534)  13.6 W12-W41 S46-S35
  

    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