Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1590 1556 (N22W56 L= 36)  1610    260        4 (1545/1441-1639/1664)  26.6 W46-W66 N16-N28
      2 1382 1534 (N27W39 L= 18)  1136    207        2 (1312/1457-1451/1595) -79.7 W35-W40 N21-N32
      3 1018 1367 (N27W08 L=348)   437    197        0 ( 967/1299-1055/1445) -90.0 W06-W09 N22-N33
      4 1181 1357 (N23W18 L=358)   253     77        2 (1149/1329-1214/1389)  11.3 W15-W20 N20-N24
      5 1357 1325 (N16W28 L=  7)   576    253        1 (1265/1253-1449/1382) -57.5 W23-W31 N10-N21
      6  884 1199 (N20E04 L=335)   191    101        0 ( 838/1182- 923/1227)   0.0 E07-E01 N19-N21
      7  620 1091 (N19E23 L=316)   496    136        2 ( 597/1022- 653/1154)  66.8 E25-E20 N15-N23
      8  826  704 (S07E19 L=321)   302    170        1 ( 796/ 643- 879/ 756)  39.8 E21-E14 S10-S05
      9 1341  696 (S20W12 L=351)   537    258        0 (1219/ 669-1420/ 750)  -8.1 W03-W17 S21-S16
     10 1529  684 (S26W24 L=  4)   831    200        2 (1453/ 628-1611/ 726)   8.1 W17-W31 S27-S23
     11  741  284 (S34E43 L=297)  2125    373        2 ( 705/ 144- 789/ 475)  32.1 E71-E28 S47-S20
  

    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