Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1102 1787 (N43E20 L=172)   990    202        1 (1060/1716-1152/1855) -33.7 E28-E13 N38-N48
      2  510 1692 (N21E55 L=137)   418     31        1 ( 487/1675- 535/1708)   0.0 E57-E52 N19-N21
      3  541 1617 (N16E48 L=144)   393     76        0 ( 528/1582- 554/1651) -45.0 E50-E45 N13-N18
      4 1240 1112 (N03W10 L=202)   271     81        1 (1197/1094-1286/1126)  11.3 W07-W12 N01-N04
      5  237 1027 (S24E54 L=139)   768    229        1 ( 172/ 974- 309/1072)  30.3 E59-E47 S27-S20
      6 1395 1039 (N03W20 L=213)   337     49        1 (1374/1006-1418/1072)  90.0 W18-W21 N01-N05
      7 1419  862 (S06W27 L=219)   975    162        3 (1395/ 773-1442/ 954) -52.1 W23-W30 S10-N00
      8 1183  801 (S15W15 L=207)  1265    293        1 (1039/ 748-1325/ 879)  42.9 W07-W21 S21-S07
      9 1849  818 (N06W62 L=254)   473    170        0 (1810/ 801-1891/ 836)  15.3 W57-W68 N04-N07
     10  866  529 (S40W05 L=197)  1690    299        4 ( 700/ 499-1043/ 563)  15.6 E08-W17 S44-S35
     11 1451  464 (S26W45 L=237)  1278    267        0 (1322/ 439-1553/ 495)  32.7 W36-W51 S29-S20
     12 1576  412 (S24W58 L=250)   370     51        0 (1558/ 386-1596/ 434) -56.3 W55-W59 S25-S22
  

    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