Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1328 1728 (N36W45 L=101)   747    109        1 (1291/1698-1368/1749)   0.0 W41-W49 N33-N37
      2 1585 1523 (N16W51 L=107)   761    208        1 (1488/1471-1662/1575) -51.3 W45-W55 N11-N21
      3 1008 1440 (N28W11 L= 67)   295     83        0 ( 972/1415-1050/1455) -11.3 W08-W13 N27-N29
      4 1122 1272 (N16W12 L= 68)   232     75        1 (1103/1234-1148/1311)  36.9 W10-W14 N13-N17
      5  156 1252 (N39E60 L=356)   865    245        0 ( 114/1233- 198/1270) -12.5 E64-E55 N37-N40
      6  853 1214 (N20E05 L= 51)   479    114        1 ( 790/1191- 912/1233) -29.7 E08-E01 N17-N21
      7 1480 1016 (S08W26 L= 82)   211     31        0 (1461/ 992-1499/1044)  33.7 W24-W27 S08-S06
      8  286  969 (N19E49 L=  7)   515    120        1 ( 243/ 941- 332/ 986)  -7.1 E53-E45 N17-N20
      9 1816  743 (S35W48 L=104)  2404    267        2 (1740/ 633-1869/ 855)  42.5 W39-W53 S40-S28
     10 1325  618 (S26W07 L= 63)  1218    331        1 (1291/ 461-1378/ 749)  60.9 W03-W13 S37-S19
     11  987  524 (S23E16 L= 40)   934    301        0 ( 909/ 419-1036/ 656)  38.7 E24-E09 S27-S15
     12  703  475 (S18E35 L= 21)   904    176        4 ( 661/ 410- 742/ 542)  56.3 E38-E30 S23-S14
  

    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