Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1439 1806 (N39W64 L=199)  1335    150        1 (1378/1766-1503/1839)   0.0 W53-W83 N37-N40
      2  841 1472 (N33E01 L=135)   759    210        2 ( 752/1451- 958/1492) -21.0 E07-W06 N30-N36
      3 1841 1387 (N04W70 L=205)   601     44        2 (1813/1363-1866/1410) -23.2 W66-W73 N02-N05
      4 1071 1386 (N23W11 L=147)   313    111        0 (1016/1361-1120/1408)   0.0 W07-W14 N21-N23
      5 1525 1349 (N10W38 L=174)   649    136        0 (1479/1286-1558/1402) -82.9 W36-W39 N05-N13
      6  181 1211 (N31E62 L= 73)  2215    220        1 ( 124/1113- 249/1293) -28.3 E79-E52 N23-N37
      7 1137 1247 (N13W11 L=147)   207     86        0 (1095/1227-1178/1272)   9.5 W08-W14 N12-N13
      8 1914 1131 (S12W69 L=204)  3414    242        7 (1874/1025-1960/1248)  63.4 W62-W79 S17-S05
      9  561 1121 (N18E26 L=109)  1753    446        0 ( 367/1064- 722/1232)   4.1 E41-E13 N15-N21
     10 1833 1146 (S09W58 L=193)   408     42        1 (1818/1115-1848/1177)  56.3 W56-W59 S10-S07
     11 1840  905 (S23W55 L=190)   762     91        0 (1815/ 863-1859/ 951)  51.3 W52-W57 S25-S20
     12 1526  778 (S23W26 L=161)  1015    250        1 (1471/ 676-1575/ 883)  45.0 W20-W30 S27-S17
     13 1221  713 (S19W05 L=140)   756    375        0 (1143/ 560-1310/ 862)  43.0 E03-W12 S27-S12
     14 1745  753 (S30W44 L=180)   646     93        0 (1717/ 720-1779/ 790)  23.2 W41-W48 S32-S28
     15  387  483 (S15E59 L= 76)  1563    247        8 ( 345/ 388- 443/ 588)  66.8 E65-E55 S22-S08
     16  849  450 (S27E26 L=110)   465     99        0 ( 831/ 398- 874/ 498)  45.0 E28-E22 S30-S24
  

    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