Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1451 1573 (N37W38 L= 31)   401     63        2 (1426/1542-1477/1603)  45.0 W36-W41 N34-N38
      2 1386 1483 (N31W29 L= 22)   582    124        2 (1325/1456-1431/1525)  -7.1 W24-W32 N29-N34
      3 1442 1365 (N23W31 L= 23)   223    106        0 (1417/1329-1477/1410) -68.2 W29-W32 N20-N25
      4 1607 1320 (N19W43 L= 35)   442    118        3 (1572/1292-1658/1355) -26.6 W40-W46 N16-N21
      5  271 1078 (N10E54 L=299)   346    101        0 ( 259/1028- 297/1114)  78.7 E54-E51 N06-N11
      6  392 1008 (N05E42 L=310)   994    198        2 ( 368/ 915- 419/1106) -90.0 E44-E40 N00-N11
      7  116 1013 (N05E75 L=277)  3113    179        1 (  93/ 924- 145/1083) -45.0 E82-E69 N00-N09
      8  867  942 (N00E10 L=342)   440    210        0 ( 808/ 870- 889/1033)  63.4 E13-E08 S03-N05
      9 1906  854 (S13W71 L= 63)   679     51        1 (1884/ 831-1932/ 877)  12.5 W67-W77 S14-S12
     10  262  747 (S11E58 L=294)  2265    350        6 ( 132/ 638- 386/ 818) -17.8 E77-E47 S17-S06
     11 1053  703 (S16W00 L=352)   197     75        1 (1030/ 668-1071/ 726) -76.0 E01-W01 S17-S13
     12  322  595 (S21E57 L=296)   930    163        1 ( 255/ 555- 396/ 636) -15.9 E63-E49 S23-S18
  

    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