Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  888 1764 (N36E29 L= 89)   770    134        1 ( 827/1727- 954/1794)  -9.5 E34-E22 N34-N37
      2 1072 1593 (N28E10 L=108)   303     74        0 (1048/1561-1109/1620)  68.2 E11-E08 N25-N30
      3  299 1578 (N14E68 L= 50)  1111    112        1 ( 267/1535- 336/1619)  71.6 E72-E65 N11-N17
      4 1412 1573 (N35W13 L=131)   809    170        1 (1323/1539-1484/1594)  14.0 W06-W18 N32-N37
      5 1126 1352 (N14E01 L=117)   217    116        0 (1094/1306-1161/1397) -38.7 E03-W02 N11-N15
      6  517 1247 (S04E35 L= 83)   357    104        0 ( 497/1197- 532/1293) -80.5 E37-E33 S06-N00
      7 1690 1163 (N16W39 L=157)   736    141        0 (1630/1106-1732/1210) -24.0 W33-W42 N13-N17
      8  414 1154 (S10E41 L= 78)   259     79        0 ( 399/1111- 426/1195)  90.0 E41-E39 S13-S07
      9  204 1155 (S13E60 L= 58)   568     97        3 ( 158/1125- 240/1179) -15.3 E66-E55 S14-S11
     10  756 1053 (S11E16 L=102)   616    233        0 ( 665/ 994- 840/1120)  48.4 E21-E12 S14-S05
     11  213  993 (S22E57 L= 61)   850    100        0 ( 167/ 956- 267/1035)  26.6 E63-E51 S24-S19
     12 1746  937 (N05W47 L=166)  1035    172        0 (1720/ 845-1763/1033) -73.3 W44-W50 N00-N10
     13  959  749 (S24W02 L=120)   310    117        0 ( 912/ 730-1009/ 787)  -9.5 E01-W05 S25-S22
     14 1751  722 (S06W55 L=173)   431     43        1 (1733/ 691-1766/ 753)  90.0 W52-W57 S08-S04
     15 1704  666 (S11W53 L=172)   846    165        3 (1647/ 604-1756/ 709)  11.3 W47-W57 S14-S07
     16 1772  628 (S11W62 L=180)   809     95        1 (1738/ 587-1813/ 669)  63.4 W58-W64 S13-S07
     17 1657  462 (S23W62 L=180)  1582    198        0 (1553/ 433-1737/ 501)  34.7 W53-W68 S27-S18
  

    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