Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  770 1821 (N37E43 L=155)   466    108        1 ( 717/1801- 811/1833)   7.1 E47-E39 N35-N38
      2  318 1501 (N07E61 L=137)   551     97        1 ( 301/1454- 341/1537) -51.3 E63-E58 N04-N09
      3 1513 1372 (N25W22 L=220)   357     95        3 (1478/1345-1551/1391)   9.5 W18-W24 N22-N26
      4 1129 1193 (N05W02 L=200)   238     75        0 (1090/1170-1169/1221)  45.0 W00-W03 N02-N07
      5  649 1121 (S10E24 L=174)   358    113        0 ( 625/1057- 664/1174) -56.3 E26-E22 S13-S06
      6  514 1015 (S18E31 L=167)   214     69        0 ( 479/ 992- 539/1043)  53.1 E33-E29 S20-S16
      7 1361  826 (S10W24 L=222)   262     92        1 (1317/ 810-1400/ 852)   0.0 W20-W26 S11-S08
      8 1469  571 (S21W40 L=238)   352    105        2 (1442/ 525-1489/ 609) -59.0 W37-W42 S23-S18
      9 1265  574 (S26W27 L=225)   227     58        1 (1242/ 558-1298/ 597)  18.4 W25-W28 S27-S25
     10  865  215 (S65W22 L=220)   444     75        1 ( 840/ 195- 889/ 242) -68.2 W19-W25 S66-S61
  

    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