Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1440 1473 (N18W39 L=326)   285    103        0 (1394/1448-1475/1503)   0.0 W35-W42 N17-N19
      2  614 1393 (N38E16 L=271)   249     72        0 ( 577/1373- 656/1407)  -8.1 E19-E12 N36-N38
      3 1263 1163 (N07W17 L=305)   833    285        3 (1139/1117-1390/1191) -21.8 W09-W24 N02-N10
      4  846 1082 (N13E08 L=279)   369    158        1 ( 768/1031- 921/1121)   0.0 E14-E03 N12-N14
      5  197  883 (N17E61 L=227)  2472    535        5 (  91/ 836- 318/ 934) -16.3 E74-E49 N11-N20
      6 1378  668 (S24W11 L=299)  1464    354        2 (1320/ 516-1445/ 815)  45.0 W03-W18 S33-S17
      7  506  621 (S05E43 L=245)   456    104        3 ( 472/ 588- 554/ 657) -51.3 E45-E39 S06-S01
      8  765  521 (S16E30 L=258)   268     43        1 ( 737/ 505- 794/ 536) -45.0 E30-E27 S17-S14
      9  656  423 (S20E41 L=247)   318    101        0 ( 614/ 409- 718/ 436) -33.7 E43-E37 S22-S18
  

    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