Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  583 1530 (N30E39 L=211)  1583    373        2 ( 521/1411- 680/1631) -28.6 E49-E27 N22-N36
      2 1591 1483 (N36W42 L=292)  1339    251        3 (1513/1369-1658/1568) -58.4 W35-W46 N28-N41
      3 1816  934 (N02W59 L=309)   759     59        3 (1786/ 906-1846/ 967)  36.9 W55-W62 N00-N03
      4 1476  929 (S01W30 L=279)   650    147        1 (1409/ 908-1555/ 963)   5.7 W25-W35 S02-N01
      5  332  787 (S19E48 L=202)   318    110        0 ( 315/ 741- 370/ 824)  50.2 E50-E44 S22-S16
      6 1812  774 (S08W62 L=312)   469     64        1 (1795/ 744-1831/ 809)  63.4 W60-W63 S10-S05
      7 1409  724 (S14W28 L=277)   960    308        1 (1298/ 657-1514/ 778) -16.4 W19-W36 S17-S10
      8 1208  682 (S18W15 L=265)   326     68        2 (1171/ 660-1235/ 707)  14.0 W12-W16 S19-S16
      9 1337  265 (S49W42 L=292)   481     26        0 (1318/ 241-1358/ 284) -35.5 W38-W45 S51-S46
  

    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