Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1610 1533 (N33W49 L= 40)   539    128        1 (1563/1503-1657/1564) -35.5 W45-W53 N30-N35
      2  757 1348 (N25E17 L=333)  3678   1107        6 ( 344/1104-1090/1551)  -7.4 E56-W04 N09-N39
      3 1440 1423 (N27W31 L= 22)   677    223        3 (1353/1394-1533/1450)   0.0 W24-W38 N25-N29
      4  358 1188 (N15E46 L=304)   702    214        1 ( 340/1101- 385/1296) -76.0 E48-E43 N10-N22
      5 1719 1225 (N13W50 L= 41)   935    228        2 (1641/1148-1791/1286) -36.9 W44-W56 N07-N16
      6  176  922 (S02E66 L=284)   764    192        0 ( 137/ 903- 222/ 939)   8.1 E73-E59 S02-N00
      7 1151  906 (S03W07 L=358)   234     57        1 (1125/ 884-1181/ 922)   0.0 W05-W09 S04-S02
      8 1479  836 (S09W29 L= 19)   997    311        3 (1376/ 747-1604/ 919) -35.2 W21-W38 S15-S03
      9  945  686 (S17E06 L=344)  1064    332        2 ( 821/ 576-1046/ 810)  43.0 E15-W00 S23-S08
     10 1123  506 (S30W05 L=355)  1592    353        2 ( 996/ 406-1282/ 594)  32.3 E04-W15 S36-S24
  

    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