Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  677 1681 (N37E30 L= 82)  3467    704       14 ( 481/1595- 887/1774) -16.3 E54-E12 N31-N45
      2  564 1347 (N14E31 L= 81)   308     86        2 ( 548/1306- 580/1389) -68.2 E32-E29 N11-N16
      3 1374 1234 (N10W21 L=133)   382    118        1 (1331/1173-1413/1292) -49.4 W17-W23 N06-N13
      4 1854 1041 (N02W59 L=171)   428     78        1 (1837/1001-1874/1078)  76.0 W57-W62 N00-N04
      5  928 1029 (S04E06 L=106)   250     68        0 ( 898/ 998- 950/1064)  63.4 E07-E04 S06-S02
      6  746  691 (S26E17 L= 95)   936    319        0 ( 630/ 582- 879/ 800)  37.9 E26-E08 S33-S18
      7 1404  675 (S23W27 L=139)   353    108        2 (1337/ 654-1451/ 697)  23.2 W22-W30 S24-S21
      8 1084  553 (S33W07 L=119)   396     98        1 (1025/ 537-1134/ 567)   7.1 W02-W10 S34-S32
      9  699  322 (S53E28 L= 84)   330     71        0 ( 681/ 295- 715/ 353)  51.3 E30-E25 S55-S50
  

    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