Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  607 1524 (N20E38 L= 43)   348     46        0 ( 583/1504- 641/1539)  33.7 E40-E36 N18-N21
      2 1416 1144 (N13W22 L=103)   661    171        2 (1345/1090-1481/1213) -20.0 W15-W26 N10-N16
      3  796  927 (S12E12 L= 70)   339    101        1 ( 782/ 878- 811/ 984) -63.4 E13-E10 S14-S08
      4 1858  913 (N09W63 L=145)   444    130        0 (1829/ 893-1891/ 926)  20.6 W59-W67 N07-N10
      5 1898  876 (N08W70 L=152)   696    140        0 (1874/ 840-1921/ 900) -11.3 W66-W76 N06-N08
      6 1565  838 (S02W38 L=120)   497    105        1 (1533/ 788-1588/ 893) -78.7 W35-W39 S04-N01
      7 1143  618 (S23W16 L= 98)   274     94        1 (1095/ 603-1193/ 633)   9.5 W13-W19 S24-S22
      8 1592  599 (S15W47 L=129)   306     45        2 (1572/ 578-1615/ 621) -45.0 W45-W49 S16-S13
      9 1646  526 (S18W56 L=137)  1737    320        3 (1580/ 435-1712/ 594)  76.0 W50-W59 S24-S12
     10 1524  549 (S19W44 L=126)   392     92        0 (1493/ 517-1568/ 586)  63.4 W42-W46 S22-S16
     11  959  569 (S30W06 L= 87)   219     94        0 ( 915/ 559-1007/ 579)  18.4 W02-W08 S31-S29
     12 1710  509 (S18W63 L=145)   555    145        0 (1687/ 483-1740/ 537) -15.9 W60-W67 S18-S16
  

    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