Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1025 1789 (N54W33 L=243)   328     55        1 ( 998/1769-1053/1807) -56.3 W31-W35 N52-N56
      2  510 1591 (N57E29 L=180)   449     51        1 ( 488/1567- 544/1613)  -7.1 E32-E24 N54-N58
      3  322 1192 (N32E47 L=163)   320     70        0 ( 287/1179- 365/1207) -20.6 E50-E42 N30-N33
      4 1399 1063 (S00W22 L=232)   209     52        0 (1374/1046-1428/1076)   0.0 W20-W24 S01-N00
      5 1598  818 (S20W30 L=240)   709    270        1 (1503/ 706-1683/ 878)  14.0 W21-W37 S24-S17
      6 1010  800 (S05E06 L=203)   200     66        0 ( 992/ 768-1029/ 834) -90.0 E06-E05 S07-S02
      7 1635  728 (S27W32 L=242)   762    210        0 (1596/ 636-1681/ 813)  48.4 W27-W36 S32-S22
      8  952  526 (S20E17 L=192)  1213    437        5 ( 808/ 414-1071/ 659)  18.4 E30-E06 S24-S13
      9 1133  474 (S28E08 L=201)   270     65        2 (1104/ 429-1183/ 509)  20.6 E11-E03 S31-S27
     10  448  416 (S17E59 L=151)   878     79        2 ( 427/ 373- 475/ 451)  45.0 E63-E54 S20-S14
     11 1320  386 (S41W02 L=211)   442     94        0 (1284/ 361-1366/ 418)   9.5 E01-W05 S42-S39
  

    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