Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1483 1543 (N33W42 L=355)  1027    214        1 (1392/1500-1580/1589) -31.6 W35-W49 N29-N37
      2 1583 1309 (N17W42 L=355)  1358    333        5 (1546/1179-1637/1427)  59.0 W36-W50 N08-N24
      3 1788 1267 (N12W60 L= 13)   534     54        2 (1759/1238-1812/1296) -53.1 W57-W62 N09-N13
      4 1660  649 (S25W44 L=357)   383    101        0 (1612/ 618-1700/ 676)  14.0 W39-W47 S26-S23
      5  790  580 (S21E19 L=294)   222     90        0 ( 748/ 557- 831/ 602)  26.6 E22-E16 S22-S19
      6  600  547 (S22E34 L=279)   294     45        1 ( 579/ 519- 623/ 574) -76.0 E35-E32 S24-S20
      7 1082  503 (S29E01 L=313)   439    176        1 (1062/ 425-1105/ 559)  78.7 E02-W01 S34-S24
      8  872  440 (S32E17 L=297)   673    232        0 ( 835/ 340- 937/ 530)  49.4 E22-E10 S39-S25
      9  771  296 (S43E30 L=283)   558     66        1 ( 733/ 268- 799/ 325) -45.0 E33-E28 S46-S40
  

    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