Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1816 1435 (N14W66 L= 67)   888    166        1 (1772/1399-1861/1476) -35.5 W61-W71 N11-N16
      2  793 1260 (N17E11 L=350)   300    135        1 ( 770/1192- 811/1326) -90.0 E11-E09 N12-N21
      3  228 1269 (N25E54 L=307)   598    179        0 ( 182/1251- 278/1286) -15.3 E59-E48 N23-N27
      4  266 1144 (N17E50 L=311)  1473    174        2 ( 224/1067- 317/1238)  60.9 E54-E45 N12-N22
      5  148  981 (N09E65 L=296)  2079    297        9 ( 101/ 911- 189/1040) -69.4 E71-E59 N04-N12
      6  772  722 (S14E20 L=342)   445    171        0 ( 739/ 646- 803/ 808)  52.1 E23-E16 S18-S09
      7  636  738 (S12E28 L=333)   280     98        0 ( 623/ 688- 647/ 786)  71.6 E29-E26 S14-S08
      8  418  652 (S14E46 L=315)   655    126        0 ( 379/ 589- 444/ 701) -81.9 E47-E43 S17-S10
      9 1269  597 (S29W10 L= 11)   873    198        2 (1234/ 520-1301/ 699)  76.0 W06-W12 S34-S22
     10  287  534 (S19E63 L=298)  1754    322        3 ( 254/ 455- 319/ 625)  54.0 E67-E56 S24-S13
     11 1049  463 (S35E07 L=354)   262     87        0 (1033/ 424-1071/ 509)  81.9 E08-E05 S38-S31
  

    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