Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  969 1762 (N49W05 L= 14)   371     47        0 ( 955/1730- 984/1795)  80.5 W03-W06 N46-N52
      2 1505 1547 (N27W40 L= 49)  1067    183        2 (1402/1518-1589/1569) -12.1 W32-W47 N24-N28
      3  813 1442 (N26E10 L=359)  1617    495        0 ( 591/1352-1009/1527) -28.2 E26-W02 N18-N34
      4 1368 1448 (N21W27 L= 36)   378     93        0 (1320/1433-1421/1462)  -8.1 W23-W30 N20-N22
      5  499 1438 (N29E33 L=336)   424     83        1 ( 454/1418- 549/1458) -20.6 E37-E29 N27-N30
      6 1806 1410 (N16W64 L= 73)   733     73        5 (1769/1381-1846/1436)   0.0 W59-W70 N14-N17
      7 1430 1216 (N07W27 L= 36)   373    130        0 (1389/1172-1481/1272) -54.5 W24-W29 N03-N10
      8 1666  970 (S10W41 L= 51)   457     95        1 (1642/ 921-1696/1021)  63.4 W39-W44 S12-S06
      9 1860  886 (S16W61 L= 70)  6252    216       18 (1764/ 756-1949/ 996)  33.7 W51-W75 S23-S09
     10 1762  652 (S30W55 L= 64)   678    148        0 (1731/ 585-1781/ 722)  49.6 W52-W56 S34-S25
     11  291  637 (S18E59 L=311)   933    165        0 ( 245/ 599- 320/ 694)  45.0 E64-E54 S19-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