Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1127 1678 (N45E04 L=120)  1296    253        4 (1086/1594-1171/1767) -66.8 E08-W00 N39-N53
      2  527 1579 (N27E46 L= 78)  2663    435        5 ( 387/1475- 630/1751) -87.3 E55-E41 N19-N40
      3 1452 1584 (N43W25 L=149)   377     94        2 (1407/1563-1495/1601)   0.0 W20-W29 N42-N44
      4 1534 1535 (N41W32 L=156)   337    116        0 (1500/1503-1569/1569)  36.9 W28-W36 N38-N44
      5  537 1417 (N17E39 L= 85)   479    142        0 ( 479/1381- 575/1473)  63.4 E42-E36 N14-N20
      6 1452 1432 (N32W24 L=148)   391    112        1 (1420/1384-1473/1469) -80.5 W21-W25 N28-N34
      7  892 1306 (N15E12 L=112)   200     91        0 ( 849/1293- 942/1321)   0.0 E15-E09 N14-N16
      8 1142 1211 (N13W05 L=129)   697    159        4 (1076/1138-1194/1275)  76.0 W00-W07 N08-N17
      9 1690 1208 (N21W43 L=167)   283     29        0 (1675/1180-1705/1230) -63.4 W41-W44 N19-N22
     10 1247 1044 (N04W13 L=137)   528    120        1 (1211/ 973-1277/1111) -54.5 W10-W15 N00-N08
     11  347  932 (S15E45 L= 79)   273     84        2 ( 333/ 875- 367/ 980)  90.0 E45-E42 S18-S11
     12 1615  920 (N02W40 L=164)   494     78        3 (1589/ 884-1635/ 963) -53.1 W37-W42 N00-N04
     13  799  896 (S11E12 L=112)   299    119        0 ( 740/ 881- 853/ 918)  23.2 E15-E08 S11-S08
     14 1155  801 (S11W11 L=135)   655    135        2 (1112/ 735-1190/ 863) -69.4 W08-W14 S15-S07
     15 1709  703 (S10W53 L=177)   549    124        0 (1672/ 661-1735/ 769)  74.1 W50-W54 S12-S05
     16  659  569 (S34E19 L=105)   241     66        0 ( 627/ 553- 682/ 589)  36.9 E21-E17 S35-S32
  

    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