Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1183 1850 (N52E19 L= 60)   410     68        3 (1144/1835-1224/1865)  38.7 E22-E15 N49-N55
      2 1098 1742 (N39E18 L= 62)   468     79        3 (1060/1716-1139/1764)  18.4 E20-E14 N37-N40
      3 1278 1671 (N39E03 L= 77)  1060    234        7 (1226/1619-1343/1720) -21.8 E07-W03 N34-N42
      4 1523 1632 (N46W18 L= 97)   450     82        1 (1490/1605-1566/1656)   0.0 W14-W22 N43-N47
      5  897 1534 (N19E21 L= 58)   719    192        0 ( 842/1464- 967/1613) -24.8 E27-E14 N15-N23
      6 1352 1428 (N24W09 L= 88)   977    359        2 (1301/1249-1396/1568) -56.3 W02-W15 N14-N33
      7  260 1487 (N04E68 L= 11)  1081    230        4 ( 232/1437- 293/1525) -31.0 E76-E63 N00-N06
      8 1772 1215 (N25W43 L=123)  1429    249        2 (1679/1133-1849/1315) -19.4 W34-W51 N22-N28
      9  246  929 (S29E51 L= 28)   650    175        2 ( 195/ 899- 303/ 955)  11.3 E56-E46 S31-S27
     10  728  453 (S49E03 L= 76)   611    136        1 ( 676/ 424- 796/ 496)   5.2 E08-W03 S51-S46
     11  864  212 (S64W26 L=105)   435     53        1 ( 838/ 193- 889/ 228)  45.0 W22-W29 S66-S62
  

    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