Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  446 1547 (N29E53 L=212)  2205    329        2 ( 401/1432- 510/1653) -33.1 E67-E42 N21-N36
      2 1416 1525 (N38W26 L=290)   620    148        2 (1366/1465-1468/1566) -41.2 W21-W29 N33-N40
      3 1128 1201 (N13W05 L=270)   268    128        1 (1069/1157-1164/1226) -31.0 W01-W07 N10-N14
      4 1668  960 (N03W44 L=309)   712    119        1 (1617/ 917-1709/1002)  45.0 W39-W47 N00-N05
      5 1207  957 (S01W12 L=276)   546    133        1 (1147/ 930-1272/ 984)  15.9 W08-W15 S02-N01
      6  245  817 (S19E58 L=207)  2015    191        4 ( 214/ 733- 284/ 896)  61.4 E62-E52 S24-S13
      7 1093  794 (S12W06 L=271)   403    103        2 (1040/ 773-1151/ 814)   0.0 W02-W09 S13-S10
      8 1277  706 (S16W19 L=284)   589    208        2 (1201/ 677-1343/ 747)   5.7 W13-W23 S17-S12
      9  960  732 (S17E01 L=263)   310    128        0 ( 910/ 710-1015/ 745)   9.5 E04-W02 S18-S16
     10  439  683 (S25E39 L=225)   631    131        1 ( 378/ 667- 494/ 707)  10.3 E45-E34 S26-S23
     11 1798  612 (S18W69 L=334)  1528    146        3 (1770/ 541-1844/ 667) -77.5 W64-W74 S23-S14
  

    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