Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  837 1647 (N51W06 L=290)   622     94        1 ( 803/1606- 877/1687)  16.7 W01-W11 N48-N53
      2 1173 1598 (N37W28 L=312)   500    107        1 (1131/1562-1217/1621) -31.0 W24-W32 N34-N38
      3 1583 1516 (N20W53 L=337)  1060    207        5 (1513/1464-1650/1551) -31.0 W46-W60 N16-N22
      4 1377 1413 (N20W32 L=317)   421    158        0 (1361/1334-1398/1480)  66.0 W30-W35 N15-N24
      5  299 1292 (N39E51 L=233)   889    146        2 ( 228/1248- 347/1333) -22.8 E63-E44 N34-N42
      6 1491 1257 (N08W34 L=318)   262     67        0 (1463/1223-1514/1290) -76.0 W33-W35 N06-N10
      7 1854 1168 (S08W62 L=347)   875    109        2 (1826/1112-1880/1212)  26.6 W57-W66 S10-S05
      8  734 1184 (N23E14 L=270)   962    365        1 ( 561/1157- 885/1204) -23.2 E25-E04 N18-N27
      9 1186 1152 (N11W12 L=297)   398    154        0 (1112/1135-1248/1171)  -6.3 W07-W16 N08-N12
     10 1504 1050 (S03W29 L=313)   766    180        3 (1461/ 968-1547/1126)  41.2 W24-W33 S07-N00
     11  268 1026 (N22E54 L=231)  1125    175        1 ( 213/ 961- 322/1084) -39.3 E60-E48 N17-N26
     12 1019 1048 (N09W00 L=285)   216    102        0 ( 975/1024-1056/1084)  11.3 E02-W03 N07-N09
     13  430 1013 (N19E39 L=246)   250     75        1 ( 397/1000- 466/1035)   0.0 E41-E35 N18-N19
     14 1522  717 (S24W24 L=308)  1371    301        1 (1443/ 599-1606/ 829)  33.7 W16-W31 S29-S18
     15  782  704 (S06E22 L=263)   211     71        0 ( 744/ 684- 814/ 723) -53.1 E23-E19 S07-S03
     16 1318  367 (S43W02 L=286)  2358    361        6 (1189/ 278-1415/ 508)  24.4 E10-W12 S48-S34
     17  792  286 (S34E39 L=245)   498     81        1 ( 754/ 262- 831/ 317)  12.5 E43-E34 S35-S32
     18 1055  197 (S50E28 L=257)   609     82        1 (1019/ 177-1100/ 216)  -8.1 E32-E23 S51-S47
  

    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