Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  798 1801 (N43E33 L=166)   432     69        1 ( 768/1776- 825/1828) -20.6 E37-E29 N40-N46
      2  711 1559 (N23E28 L=171)   247     74        1 ( 688/1531- 738/1592) -45.0 E30-E25 N21-N25
      3  316 1523 (N18E61 L=138)  2663    158        7 ( 249/1466- 377/1589)  28.6 E68-E55 N14-N22
      4 1075 1503 (N24E03 L=196)   567    103        1 (1035/1447-1098/1556) -81.9 E05-E01 N20-N27
      5  308 1401 (N10E56 L=144)  1462    177        4 ( 255/1322- 361/1457)  60.9 E60-E51 N04-N13
      6 1752 1290 (N22W47 L=247)  2354    413        2 (1587/1147-1847/1402) -26.6 W33-W57 N14-N26
      7 1147 1356 (N15W03 L=202)   456     99        2 (1092/1321-1195/1384)   0.0 W00-W06 N13-N17
      8  516 1239 (N01E34 L=165)   505    148        0 ( 472/1211- 562/1276)  26.6 E37-E31 N00-N03
      9  911 1203 (N03E09 L=190)   364    130        0 ( 871/1150- 948/1245) -38.7 E11-E06 N00-N05
     10 1353  914 (S08W21 L=220)   264     68        0 (1332/ 891-1380/ 950) -45.0 W19-W23 S08-S05
     11  661  805 (S23E21 L=179)   542    137        1 ( 599/ 775- 727/ 834)   0.0 E24-E15 S24-S21
     12 1286  751 (S18W20 L=219)   735    299        1 (1148/ 703-1416/ 794)  28.1 W12-W27 S22-S14
     13 1687  418 (S31W70 L=269)  3070    243        3 (1587/ 365-1788/ 501)  42.7 W61-W77 S36-S23
  

    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