Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  882 1686 (N48E15 L=123)  2444    425        5 ( 787/1557-1024/1743)  32.9 E22-E01 N37-N53
      2 1364 1634 (N44W28 L=166)  1514    354        2 (1272/1560-1487/1682)  22.2 W19-W41 N38-N48
      3 1577 1605 (N42W47 L=185)   438     56        1 (1554/1583-1605/1622) -14.0 W44-W49 N40-N43
      4 1654 1545 (N37W52 L=190)   438     81        0 (1626/1528-1681/1562)   9.5 W48-W54 N36-N38
      5  446 1421 (N26E44 L= 95)   437    191        2 ( 434/1348- 470/1496) -63.4 E45-E40 N21-N31
      6 1457 1458 (N31W31 L=169)   284     62        0 (1430/1432-1485/1489) -36.9 W29-W33 N29-N33
      7 1309 1429 (N29W19 L=158)   230     72        0 (1290/1400-1329/1468) -76.0 W18-W20 N27-N31
      8  991 1304 (N21E03 L=135)   424     90        1 ( 950/1270-1049/1334) -26.6 E05-W01 N18-N22
      9  482 1214 (N13E37 L=101)   264    115        0 ( 458/1160- 506/1269) -50.2 E39-E34 N10-N16
     10  945 1223 (N15E05 L=133)   488    212        3 ( 839/1202-1049/1248)   8.7 E12-W01 N13-N17
     11  830 1168 (N12E12 L=126)   222     52        2 ( 813/1140- 847/1199) -76.0 E13-E11 N09-N13
     12 1584 1067 (N06W37 L=175)   279     86        0 (1554/1037-1618/1103) -45.0 W34-W39 N04-N08
     13 1281  905 (S04W16 L=154)   299    116        0 (1266/ 846-1308/ 951)  74.1 W15-W17 S07-N00
     14  288  822 (S12E53 L= 86)  1005    298        3 ( 226/ 746- 340/ 899) -38.7 E59-E47 S16-S07
     15  992  688 (S18E01 L=137)   515    159        2 ( 954/ 622-1025/ 763)  83.7 E03-W00 S22-S13
     16  485  396 (S42E45 L= 93)  1311    257        1 ( 433/ 356- 528/ 460)  33.7 E53-E38 S45-S35
     17 1660  446 (S35W54 L=193)   492    172        1 (1619/ 434-1709/ 457)  -5.7 W49-W59 S36-S34
     18 1314  376 (S40W25 L=163)  1539    344        2 (1174/ 329-1485/ 429)  11.8 W14-W39 S44-S36
  

    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