Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  501 1496 (N31E43 L= 37)   673    117        1 ( 446/1476- 563/1513)  -4.8 E48-E36 N29-N32
      2  904 1363 (N24E09 L= 71)   475    165        0 ( 834/1352- 999/1377)   5.2 E13-E02 N23-N24
      3 1479 1199 (N15W30 L=109)   315    115        0 (1420/1181-1517/1219)  -8.1 W25-W32 N13-N15
      4  147  895 (S09E71 L=  8)  1637    134        2 ( 116/ 846- 175/ 954) -54.5 E78-E67 S12-S05
      5 1401  840 (S08W25 L=104)   260     55        1 (1375/ 812-1427/ 862)  45.0 W22-W26 S09-S06
      6 1648  735 (S14W45 L=124)  1031    194        2 (1587/ 671-1712/ 816)  41.6 W40-W49 S18-S09
      7  365  704 (S20E47 L= 32)   969    188        1 ( 334/ 620- 414/ 785) -79.7 E49-E44 S25-S14
      8 1056  496 (S31W04 L= 84)   581    144        0 (1038/ 428-1081/ 560) -73.3 W02-W06 S36-S26
      9 1105  375 (S41W09 L= 89)   362     57        0 (1080/ 350-1128/ 403)  53.1 W06-W10 S42-S38
  

    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