Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1485 1710 (N57W19 L=216)  1100    119        4 (1436/1677-1531/1740)   0.0 W13-W25 N53-N60
      2  562 1668 (N25E52 L=145)   720    104        2 ( 530/1628- 594/1714) -60.3 E55-E49 N21-N28
      3 1265 1631 (N41E00 L=196)   469     71        1 (1232/1590-1296/1663) -21.8 E03-W02 N38-N43
      4 1370 1464 (N32W11 L=208)   248     56        1 (1347/1438-1393/1489) -33.7 W09-W13 N30-N33
      5  148 1111 (S16E70 L=127)   910    108        3 ( 122/1059- 173/1151) -29.1 E76-E65 S19-S13
      6  451  844 (S27E34 L=163)   782    186        1 ( 380/ 783- 491/ 904)  51.3 E39-E30 S32-S22
      7 1660  722 (S04W49 L=246)   333     45        1 (1635/ 703-1688/ 745)  26.6 W47-W51 S05-S03
      8  446  679 (S38E33 L=164)   671     96        2 ( 416/ 643- 499/ 715)  29.7 E36-E28 S40-S35
      9 1577  629 (S12W46 L=243)   404     63        1 (1537/ 609-1609/ 656)  31.0 W43-W48 S13-S09
     10  501  517 (S48E28 L=169)   662    103        1 ( 469/ 474- 539/ 558)  54.5 E31-E23 S51-S44
     11  598  357 (S59E17 L=179)   701     66        1 ( 569/ 324- 630/ 393)  56.3 E21-E13 S63-S54
     12 1573  345 (S28W68 L=265)  1101    145        2 (1502/ 321-1633/ 375)  21.0 W60-W75 S30-S25
  

    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