Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1338 1735 (N47W42 L=131)  1995    304        4 (1275/1644-1391/1811)  52.4 W35-W50 N39-N52
      2  920 1555 (N41W00 L= 89)   359     95        2 ( 903/1519- 935/1592) -90.0 E01-W01 N38-N44
      3 1100 1500 (N35W13 L=102)   903    176        1 (1026/1428-1162/1557) -65.6 W07-W15 N28-N39
      4 1470 1501 (N28W39 L=129)   280     61        1 (1443/1471-1489/1525) -76.0 W38-W41 N26-N30
      5  857 1431 (N33E06 L= 83)   414     76        1 ( 836/1389- 881/1476) -71.6 E07-E04 N30-N36
      6 1692 1246 (N09W49 L=138)  1258    297        4 (1565/1162-1804/1326) -37.9 W40-W58 N02-N16
      7 1265 1200 (N13W17 L=106)   295    101        0 (1212/1169-1310/1224)   8.1 W13-W20 N12-N14
      8  347 1040 (N13E46 L= 43)   558    147        0 ( 299/1002- 394/1095) -36.9 E50-E42 N10-N17
      9 1548  575 (S31W31 L=120)  2024    392        2 (1373/ 490-1689/ 698)  19.8 W17-W42 S34-S24
     10  766  488 (S25E25 L= 64)   824    108        3 ( 700/ 458- 834/ 533)  20.0 E30-E19 S26-S22
  

    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