Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  568 1733 (N34E53 L= 89)   455     88        1 ( 543/1713- 594/1749)  56.3 E56-E49 N32-N35
      2 1277 1582 (N39W06 L=148)   267     52        1 (1254/1559-1305/1599)  45.0 W03-W07 N37-N40
      3 1385 1529 (N38W15 L=157)   433    104        1 (1332/1502-1436/1551)   6.3 W10-W19 N35-N39
      4  917 1499 (N25E17 L=126)   378    135        1 ( 865/1451- 972/1550) -24.0 E21-E12 N23-N28
      5 1628 1407 (N35W36 L=178)  1307    235        1 (1549/1353-1723/1470) -14.9 W28-W44 N32-N37
      6 1761 1337 (N33W48 L=191)   986    188        0 (1703/1308-1807/1362)  20.0 W42-W53 N30-N35
      7 1413 1173 (N15W21 L=164)  1011    325        0 (1369/1042-1491/1253) -81.3 W17-W25 N07-N20
      8 1164  975 (S02W09 L=152)   426    140        1 (1082/ 947-1211/1021)  39.8 W04-W11 S04-N01
      9  512  700 (S31E29 L=114)   309     58        0 ( 490/ 669- 538/ 728)  68.2 E30-E26 S33-S28
     10 1166  689 (S18W15 L=158)   238     37        0 (1143/ 665-1188/ 710) -26.6 W13-W17 S19-S17
     11 1629  655 (S11W48 L=191)   703    161        2 (1575/ 628-1686/ 690)  -5.7 W43-W53 S12-S09
     12 1255  250 (S47W42 L=185)   369     38        0 (1231/ 233-1274/ 269)  45.0 W40-W44 S47-S44
  

    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