Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  499 1402 (N33E36 L=319)   791    181        1 ( 424/1337- 570/1444) -28.1 E44-E29 N28-N36
      2 1147 1241 (N17W10 L=  6)   191     57        0 (1134/1210-1161/1272) -76.0 W09-W10 N15-N19
      3  173 1227 (N23E70 L=286)   585     72        1 ( 158/1192- 189/1253)  63.4 E73-E67 N20-N24
      4 1725 1054 (S01W48 L= 44)   900    208        0 (1702/ 970-1754/1143)  66.0 W45-W52 S06-N03
      5 1677  850 (S13W43 L= 38)   322     76        0 (1644/ 829-1702/ 879)  21.8 W39-W44 S13-S11
      6 1836  817 (S18W60 L= 55)   766    133        0 (1803/ 780-1863/ 859)  59.0 W55-W63 S20-S15
      7  629  749 (S08E28 L=328)   378     64        1 ( 599/ 725- 664/ 772)  14.0 E29-E25 S09-S06
      8  548  682 (S12E35 L=321)   405    104        0 ( 525/ 639- 573/ 720)  59.0 E36-E32 S14-S09
      9 1257  669 (S19W11 L=  7)   316    125        1 (1230/ 631-1286/ 715) -90.0 W09-W13 S21-S16
     10  732  647 (S15E23 L=333)   490    119        0 ( 687/ 618- 762/ 700)  45.0 E25-E20 S17-S11
     11 1428  568 (S28W23 L= 19)  1695    452        2 (1284/ 461-1616/ 691)  23.7 W12-W38 S34-S22
     12  934  640 (S18E10 L=346)   485    155        0 ( 854/ 624-1004/ 654)   0.0 E14-E05 S18-S16
     13 1737  629 (S29W52 L= 47)   554     72        1 (1699/ 612-1775/ 648)  -6.3 W47-W56 S30-S27
     14  993  550 (S24E07 L=348)   429    111        0 ( 968/ 502-1034/ 598)  45.0 E09-E03 S27-S21
     15 1485  469 (S37W29 L= 25)  1089    263        0 (1367/ 430-1586/ 531)  12.5 W19-W38 S38-S33
     16 1017  370 (S38E08 L=347)   431     87        1 ( 977/ 352-1074/ 388)   0.0 E11-E03 S39-S36
     17 1078  303 (S45E05 L=350)  1816    243        4 (1002/ 224-1167/ 374)  42.9 E12-W03 S52-S39
  

    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