Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  939 1468 (N35W05 L=298)  1198    338        1 ( 794/1442-1094/1498) -13.4 E05-W16 N31-N38
      2 1292 1427 (N24W27 L=320)   773    272        2 (1172/1341-1381/1471) -50.7 W20-W30 N17-N29
      3 1276 1265 (N15W21 L=313)   619    204        1 (1234/1203-1333/1309)  71.6 W17-W23 N11-N17
      4 1473 1269 (N11W33 L=326)   220     78        0 (1438/1252-1508/1286) -36.9 W31-W35 N09-N12
      5 1767 1055 (S09W50 L=342)   646    238        1 (1702/ 972-1811/1160)  33.7 W42-W55 S12-S04
      6  649  946 (N10E24 L=269)  1009    344        1 ( 614/ 799- 692/1099)  87.0 E26-E21 N00-N19
      7 1913 1052 (S15W68 L=  1)   777    113        4 (1887/1005-1941/1093)  53.1 W63-W75 S17-S12
      8  226 1011 (N20E59 L=234)   701    128        0 ( 195/ 964- 276/1061) -32.0 E62-E53 N16-N23
      9  149  971 (N18E70 L=223)   815     70        2 ( 131/ 927- 175/1004)  63.4 E73-E65 N15-N19
     10  280  591 (S06E65 L=228)  6833    503       10 ( 198/ 423- 376/ 778)  60.5 E81-E50 S18-N05
     11 1678  751 (S25W38 L=330)   362     74        1 (1652/ 725-1712/ 774)  11.3 W35-W40 S26-S23
     12 1182  643 (S19W01 L=294)   259     56        2 (1145/ 622-1207/ 663)   0.0 E01-W02 S20-S17
  

    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