Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1411 1554 (N31W39 L= 19)   497     55        1 (1379/1527-1450/1573)   9.5 W36-W42 N29-N32
      2  693 1337 (N31E17 L=322)  1774    468        4 ( 497/1258- 881/1444) -32.7 E33-E05 N22-N40
      3 1147 1137 (N11W09 L=349)   282     64        0 (1125/1102-1170/1173)  53.1 W07-W11 N09-N13
      4  466 1046 (N15E36 L=304)   301    111        0 ( 434/1008- 517/1074)  23.2 E38-E31 N13-N16
      5  904 1035 (N09E07 L=333)   195     52        0 ( 894/1005- 919/1068)  76.0 E07-E06 N07-N11
      6  247 1020 (N16E56 L=284)   632     88        1 ( 208/ 993- 292/1055) -18.4 E60-E51 N13-N18
      7 1209  807 (S09W08 L=347)   661    137        0 (1189/ 738-1230/ 883)  71.6 W05-W09 S13-S04
      8 1439  728 (S18W21 L=  1)   314    106        0 (1391/ 709-1485/ 752)  -8.1 W17-W24 S19-S16
      9  678  599 (S14E29 L=311)   828    279        1 ( 628/ 543- 743/ 663) -76.0 E31-E24 S17-S09
     10  583  557 (S16E37 L=303)   651    186        2 ( 516/ 490- 619/ 601) -76.0 E40-E34 S20-S12
     11  898  395 (S31E21 L=319)  1135    233        1 ( 813/ 338- 986/ 485)  25.2 E28-E11 S34-S25
     12  723  292 (S37E39 L=300)   332     38        1 ( 704/ 274- 747/ 315)  11.3 E42-E37 S38-S34
  

    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