Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1699 1607 (N22W73 L= 18)  1285     80        4 (1672/1574-1738/1643) -45.0 W68-W79 N19-N24
      2  769 1415 (N36E08 L=298)  1010    255        0 ( 644/1387- 881/1437) -19.4 E17-W00 N31-N39
      3 1112 1370 (N25W14 L=319)   653    266        0 (1003/1293-1201/1401) -48.0 W07-W17 N18-N29
      4 1068 1192 (N16W06 L=312)   319    137        1 (1048/1134-1101/1240)  63.4 W04-W09 N12-N18
      5 1827 1031 (S12W55 L=  1)   475     84        1 (1809/ 992-1847/1074)  78.7 W53-W57 S14-S09
      6  857 1034 (N11E10 L=296)   216    102        0 ( 812/1000- 897/1061) -51.3 E11-E07 N08-N13
      7  457 1020 (N17E36 L=269)   326     52        2 ( 429/ 996- 492/1044) -21.8 E38-E33 N15-N18
      8  230  895 (N13E60 L=246)  1726    154        5 ( 192/ 830- 271/ 951)  41.2 E65-E54 N08-N16
      9 1270  875 (S07W11 L=317)   481    216        3 (1222/ 807-1344/ 928)  35.5 W07-W15 S09-S04
     10  458  884 (N09E39 L=267)   448     73        1 ( 437/ 841- 488/ 922) -90.0 E40-E36 N06-N11
     11  205  700 (N01E69 L=236)   748     83        2 ( 172/ 673- 240/ 736) -53.1 E75-E64 S01-N03
     12  521  374 (S24E54 L=252)   559     58        1 ( 488/ 358- 554/ 394) -18.4 E57-E50 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