Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1320 1264 (N13W12 L=144)  1090    296        1 (1281/1111-1354/1407) -58.0 W07-W17 N04-N20
      2  288 1325 (S01E56 L= 75)  1230    231        2 ( 216/1222- 360/1398)  48.0 E60-E50 S07-N03
      3  315 1148 (S11E49 L= 83)   466    102        0 ( 295/1095- 347/1196)  71.6 E50-E47 S14-S08
      4 1572 1091 (N09W32 L=163)   285    117        0 (1539/1035-1591/1151)  82.9 W30-W32 N04-N12
      5 1844 1013 (N12W56 L=187)  1873    372        7 (1811/ 886-1879/1112) -81.9 W52-W59 N04-N18
      6 1769 1019 (N10W48 L=180)   334     40        0 (1749/ 993-1785/1043) -71.6 W46-W50 N08-N11
      7 1389  570 (S25W34 L=166)   945    250        1 (1279/ 545-1510/ 593)  25.0 W26-W41 S27-S20
      8 1123  550 (S33W18 L=149)   204     46        0 (1094/ 539-1153/ 563)  33.7 W15-W19 S33-S31
      9 1700  487 (S21W62 L=194)  1268    246        2 (1624/ 448-1786/ 554)  42.3 W56-W69 S24-S14
  

    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