Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  984 1813 (N52E25 L=200)   446     57        0 ( 963/1791-1014/1833)   0.0 E27-E21 N49-N54
      2 1050 1722 (N46E14 L=211)  1149    146        4 ( 993/1682-1110/1773)  11.3 E17-E07 N42-N50
      3 1061 1605 (N37E09 L=216)   555     82        1 (1020/1574-1102/1650)  81.9 E10-E05 N33-N40
      4  372 1541 (N20E60 L=165)  1317    245        1 ( 302/1446- 446/1620) -15.3 E80-E47 N15-N24
      5  283 1381 (N09E60 L=164)  1732    222        0 ( 237/1293- 323/1482) -76.0 E64-E56 N03-N15
      6 1352 1172 (N13W18 L=242)   812    181        1 (1283/1126-1410/1214)  45.0 W13-W21 N09-N16
      7 1301 1103 (N08W15 L=240)   967    266        3 (1226/1028-1395/1208) -22.6 W09-W21 N04-N13
      8 1075 1115 (N06W02 L=226)   467    135        1 (1057/1050-1095/1186) -66.8 W00-W03 N02-N09
      9  407 1029 (S10E40 L=185)   715    228        1 ( 393/ 910- 429/1119) -80.5 E41-E38 S16-S04
     10 1358 1010 (N04W20 L=245)   197     69        1 (1330/ 982-1379/1044)  76.0 W18-W21 N01-N05
     11 1571 1007 (N07W35 L=259)   313    130        0 (1522/ 966-1607/1033)  -9.5 W31-W37 N05-N08
     12  319  853 (S22E47 L=178)  1062    192        1 ( 227/ 813- 404/ 879)  14.7 E57-E38 S25-S19
     13 1285  708 (S15W21 L=246)   517    160        0 (1248/ 643-1325/ 773)  77.5 W19-W22 S19-S10
     14 1677  606 (S15W54 L=278)   344     63        0 (1638/ 586-1702/ 625)  -7.1 W49-W57 S15-S13
     15 1559  436 (S27W53 L=277)   843    177        0 (1474/ 397-1634/ 474)  26.6 W46-W58 S31-S23
     16  497  418 (S50E36 L=189)   606     87        1 ( 448/ 404- 543/ 438)  13.0 E43-E30 S52-S48
     17  601  388 (S51E24 L=200)   451     44        0 ( 580/ 362- 626/ 420) -78.7 E26-E21 S52-S47
     18  698  308 (S56E14 L=210)   920    112        1 ( 663/ 275- 746/ 353) -54.5 E18-E08 S58-S51
     19  970  236 (S56W17 L=241)  1104    152        2 ( 928/ 186-1026/ 303) -27.9 W08-W26 S60-S50
  

    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