Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  394 1615 (N20E62 L=336)   597     30        0 ( 368/1597- 418/1635)  36.9 E65-E59 N18-N21
      2  470 1490 (N13E48 L=350)  1679    405        3 ( 413/1380- 569/1612) -31.0 E57-E36 N08-N20
      3 1025 1402 (N20E09 L= 30)   675    177        1 ( 973/1328-1083/1472) -45.0 E13-E03 N15-N23
      4 1447 1333 (N26W20 L= 58)   397    111        0 (1412/1275-1479/1380) -39.8 W16-W22 N22-N27
      5 1054 1212 (N09E02 L= 36)   298     56        1 (1025/1182-1078/1239) -45.0 E04-W00 N07-N10
      6  425 1194 (S05E42 L=356)  1096    189        3 ( 349/1141- 507/1235)  38.7 E46-E36 S09-N00
      7 1772 1078 (N18W49 L= 87)  1398    339        1 (1724/ 938-1805/1184) -63.4 W43-W54 N10-N24
      8  667 1026 (S10E21 L= 17)   753    233        0 ( 559/ 985- 769/1064)  15.9 E28-E14 S11-S07
      9 1150  593 (S25W18 L= 56)   382     98        0 (1093/ 577-1201/ 607)  15.9 W14-W21 S25-S23
     10  670  546 (S40E13 L= 25)   415    113        0 ( 612/ 533- 722/ 562)  26.6 E17-E09 S42-S37
  

    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