Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  511 1716 (N22E57 L= 67)   946    209        3 ( 476/1656- 534/1769) -71.6 E61-E54 N17-N26
      2 1110 1721 (N38E17 L=107)   796    115        2 (1048/1689-1170/1746)  12.5 E20-E11 N35-N39
      3  491 1628 (N15E53 L= 71)  1049    211        3 ( 400/1579- 576/1685)  45.0 E59-E48 N11-N20
      4 1222 1617 (N34E05 L=118)   725    144        3 (1173/1559-1267/1666) -40.6 E09-E01 N30-N36
      5  627 1593 (N16E41 L= 83)   384     96        1 ( 599/1560- 670/1618)  -8.1 E43-E36 N14-N17
      6 1580 1107 (N13W30 L=154)   240     63        0 (1550/1088-1609/1123)  14.0 W28-W32 N12-N13
      7 1828  827 (N06W59 L=182)   428     66        1 (1805/ 795-1844/ 863)  90.0 W57-W60 N03-N08
      8 1151  694 (S22W17 L=140)   437    161        0 (1064/ 675-1232/ 712)  21.8 W11-W21 S24-S19
      9 1639  494 (S17W58 L=181)   888    144        1 (1588/ 435-1699/ 550)  77.5 W55-W60 S21-S12
     10 1692  465 (S16W65 L=189)  1674    229        4 (1638/ 413-1748/ 508)  53.1 W60-W70 S20-S12
     11 1375  481 (S27W39 L=163)   473    109        2 (1316/ 457-1416/ 503)  14.0 W34-W42 S27-S24
     12 1480  412 (S27W50 L=174)   303     67        0 (1447/ 386-1509/ 436)  76.0 W49-W51 S28-S24
  

    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