Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  694 1742 (N61E01 L=123)   324     56        0 ( 675/1712- 716/1764)  45.0 E03-W02 N58-N62
      2 1131 1568 (N33W23 L=147)   594    120        0 (1068/1535-1181/1593) -50.2 W19-W26 N29-N35
      3 1354 1503 (N22W34 L=158)   408    128        1 (1304/1463-1407/1537)   5.7 W29-W39 N21-N23
      4  607 1184 (N24E21 L=103)   265    113        0 ( 587/1134- 637/1231) -80.5 E21-E20 N20-N26
      5  206 1125 (N29E58 L= 66)  1186    112        3 ( 169/1081- 249/1169)  39.8 E63-E52 N26-N31
      6 1837  953 (S22W54 L=177)   779     95        1 (1816/ 910-1860/ 997) -90.0 W51-W57 S25-S19
      7 1682  650 (S36W36 L=159)  2650    628        4 (1457/ 508-1892/ 856)   8.5 W14-W61 S40-S30
      8  748  670 (S08E25 L= 99)   513    133        1 ( 722/ 618- 784/ 731)  45.0 E27-E22 S11-S05
      9 1150  402 (S36E10 L=113)   239     64        0 (1119/ 382-1174/ 426)  21.8 E13-E08 S36-S34
     10  806  389 (S27E33 L= 91)   808    131        2 ( 752/ 370- 872/ 417) -35.5 E36-E28 S29-S24
     11  633  352 (S25E47 L= 77)   333     76        0 ( 595/ 340- 668/ 366)   0.0 E51-E43 S26-S24
     12  994  144 (S53E43 L= 81)   454     28        0 ( 974/ 123-1012/ 166)  50.2 E47-E38 S55-S49
  

    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