Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1515 1735 (N59W46 L=181)   741     76        3 (1479/1709-1559/1751)   9.0 W38-W57 N55-N60
      2 1405 1731 (N56W28 L=162)   392     53        0 (1385/1704-1431/1748)  39.8 W24-W32 N53-N58
      3  892 1642 (N39E20 L=115)   421     91        2 ( 850/1606- 925/1666) -23.2 E24-E17 N36-N40
      4 1400 1626 (N46W23 L=158)   321     40        0 (1374/1609-1425/1641) -11.3 W20-W25 N44-N47
      5  626 1451 (N22E33 L=102)   617    149        1 ( 586/1397- 652/1527) -77.5 E35-E31 N18-N27
      6 1181 1461 (N30W05 L=140)   270     95        0 (1165/1428-1200/1509) -59.0 W03-W06 N28-N33
      7  448 1438 (N19E47 L= 88)   628     99        4 ( 411/1401- 494/1479)  39.8 E51-E43 N16-N21
      8  550 1416 (N19E38 L= 97)   325    102        1 ( 530/1372- 574/1456) -78.7 E39-E36 N16-N21
      9 1716 1419 (N34W52 L=187)   924    150        1 (1653/1374-1779/1450) -17.1 W45-W58 N31-N35
     10  232 1191 (N01E60 L= 74)   377     71        1 ( 206/1165- 266/1208)  18.4 E63-E56 N00-N02
     11  427  821 (S19E38 L= 96)   912    149        4 ( 387/ 755- 485/ 864)  41.2 E42-E33 S23-S15
     12 1113  428 (S37W16 L=150)   559    149        1 (1053/ 381-1179/ 479)  48.8 W11-W19 S40-S32
  

    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