Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  643 1391 (N04E32 L=179)   339    120        3 ( 608/1367- 703/1413)  18.4 E34-E28 N02-N05
      2 1378 1332 (N20W12 L=222)   494     91        3 (1343/1295-1412/1363) -31.0 W09-W14 N18-N22
      3  415 1203 (S11E42 L=168)   281     53        1 ( 392/1175- 440/1225)  76.0 E44-E40 S13-S09
      4 1112  993 (S06W06 L=216)   480    163        0 (1022/ 968-1188/1009)  16.7 W00-W10 S08-S04
      5  920  928 (S15E03 L=207)   558    246        0 ( 854/ 819- 966/1003) -59.0 E08-W00 S21-S11
      6 1312  825 (S10W22 L=232)   338    111        0 (1257/ 803-1359/ 849)  33.7 W18-W24 S12-S07
      7  501  778 (S34E27 L=183)   364     80        0 ( 457/ 754- 535/ 802)   0.0 E30-E24 S35-S32
      8  819  711 (S30E03 L=207)   216     92        0 ( 774/ 682- 865/ 745)  -8.1 E07-W00 S30-S28
      9 1411  563 (S21W38 L=248)  2617    771        6 (1087/ 521-1711/ 620)  19.2 W16-W62 S28-S10
     10 1665  419 (S19W66 L=276)  1305    265        1 (1626/ 360-1684/ 467) -45.0 W59-W71 S21-S16
     11 1493  407 (S26W52 L=262)  1148    327        2 (1346/ 383-1608/ 429)  24.2 W40-W60 S30-S20
     12 1405  274 (S36W55 L=266)   855    190        1 (1344/ 238-1454/ 302)   9.5 W49-W61 S37-S34
  

    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