Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  480 1579 (N16E49 L=227)   300     43        0 ( 455/1561- 504/1596)  33.7 E50-E46 N14-N16
      2 1166 1466 (N22E01 L=275)   398     79        1 (1128/1422-1196/1487)  90.0 E03-W00 N18-N23
      3 1237 1379 (N18W05 L=281)   587    147        3 (1198/1336-1289/1459) -32.0 W00-W08 N15-N22
      4 1215 1186 (N06W07 L=284)   359     99        1 (1170/1148-1251/1223) -18.4 W04-W10 N04-N07
      5  181 1168 (S13E63 L=214)   404    110        0 ( 153/1135- 211/1208)  38.7 E65-E60 S14-S10
      6 1749  981 (N08W47 L=323)   667    157        0 (1718/ 920-1802/1018)  80.5 W44-W51 N04-N10
      7  246  896 (S29E52 L=225)   374    103        0 ( 215/ 868- 277/ 921)  26.6 E54-E48 S30-S26
      8  305  818 (S33E46 L=230)   434    113        0 ( 257/ 805- 346/ 830)  14.0 E50-E42 S33-S31
      9 1723  542 (S16W62 L=338)   488     51        0 (1695/ 509-1746/ 574)  68.2 W59-W63 S18-S13
     10 1476  472 (S27W45 L=321)  1142    245        0 (1363/ 441-1586/ 511)  34.7 W38-W51 S31-S22
     11  861  210 (S66W19 L=296)   539     31        1 ( 837/ 186- 887/ 232)   9.5 W15-W23 S68-S63
  

    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