Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  502 1664 (N58E45 L= 28)   584     82        1 ( 476/1634- 537/1691)  51.3 E50-E41 N55-N60
      2 1718 1458 (N19W62 L=135)   434     86        0 (1678/1428-1748/1486) -51.3 W59-W64 N16-N21
      3  713 1439 (N37E16 L= 58)   242     57        1 ( 679/1421- 745/1451) -18.4 E18-E12 N35-N38
      4  778 1187 (N20E13 L= 60)   363    121        0 ( 737/1162- 823/1211) -31.0 E15-E10 N17-N21
      5  396 1064 (N17E42 L= 31)   447     90        1 ( 367/1020- 420/1100) -59.0 E44-E40 N14-N19
      6 1509  721 (S21W27 L=100)   710    179        0 (1456/ 650-1564/ 790)  37.9 W22-W31 S24-S16
      7 1214  624 (S21W06 L= 79)   827    321        1 (1111/ 542-1313/ 741)  31.0 E02-W13 S25-S15
      8  713  586 (S16E27 L= 46)   445    105        1 ( 690/ 533- 749/ 628)  50.2 E29-E24 S19-S13
      9 1262  510 (S30W07 L= 80)   784    309        1 (1140/ 447-1372/ 587)  19.4 E02-W15 S32-S26
     10 1672  545 (S37W43 L=116)   476     83        0 (1636/ 525-1711/ 574)   7.1 W39-W47 S38-S35
     11 1448  459 (S38W22 L= 95)   619    211        0 (1373/ 402-1532/ 518)  24.8 W15-W28 S41-S35
     12  813  446 (S27E25 L= 48)   605    121        1 ( 786/ 403- 851/ 493) -80.5 E26-E22 S30-S23
  

    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