Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1389 1707 (N43W48 L= 13)   505     76        0 (1360/1683-1422/1727)   0.0 W43-W52 N41-N44
      2  605 1437 (N38E25 L=299)   455    132        2 ( 541/1404- 659/1463) -15.3 E31-E20 N34-N39
      3 1548 1288 (N14W39 L=  3)   837    186        1 (1515/1210-1590/1371)  84.3 W36-W40 N08-N18
      4 1141 1229 (N17W10 L=335)   371    155        1 (1067/1197-1208/1259) -35.5 W06-W14 N14-N19
      5  409 1163 (N22E41 L=284)   400    122        1 ( 370/1118- 460/1203) -36.9 E44-E36 N18-N24
      6  820 1150 (N16E11 L=314)   279     56        0 ( 800/1121- 838/1181) -90.0 E12-E10 N14-N18
      7 1715  961 (S09W46 L= 10)   312     85        1 (1695/ 926-1745/ 984)  26.6 W43-W48 S11-S08
      8 1777  946 (S11W51 L= 16)   355     93        0 (1748/ 910-1794/ 983)  36.9 W48-W53 S13-S09
      9 1561  751 (S19W31 L=356)   262     78        1 (1549/ 710-1579/ 791)  68.2 W30-W33 S21-S16
     10 1687  703 (S25W43 L=  7)   799    193        1 (1634/ 629-1739/ 777)  41.6 W38-W47 S29-S21
     11  490  399 (S28E53 L=271)  2241    249        3 ( 459/ 288- 534/ 496)  28.0 E75-E43 S38-S21
     12  960  323 (S40E18 L=307)   317     54        0 ( 935/ 304- 992/ 339)  26.6 E20-E14 S41-S38
  

    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