Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1490 1717 (N54W12 L=249)   513     71        1 (1462/1688-1513/1738) -45.0 W09-W15 N50-N56
      2  528 1650 (N17E52 L=185)   481    153        0 ( 495/1603- 548/1691) -29.1 E56-E47 N14-N19
      3  469 1628 (N15E55 L=182)   435     65        1 ( 442/1604- 501/1646)  45.0 E57-E52 N13-N16
      4  404 1524 (N07E55 L=182)  1361    236        3 ( 328/1488- 496/1574)  30.3 E62-E48 N04-N11
      5 1004 1512 (N20E15 L=222)   339     82        2 ( 969/1484-1050/1536)  21.8 E16-E11 N17-N21
      6 1323 1359 (N20W08 L=245)   349    134        0 (1291/1288-1343/1421) -45.0 W04-W10 N16-N22
      7  206 1251 (S12E63 L=174)  1920    218        3 ( 152/1162- 240/1356)  85.2 E66-E58 S17-S05
      8  717 1180 (S06E21 L=216)   417    127        2 ( 648/1156- 771/1196)  23.2 E25-E17 S07-S04
      9  492  831 (S31E29 L=208)   848    184        1 ( 424/ 806- 578/ 878)   4.8 E34-E22 S32-S29
     10  852  752 (S27E03 L=234)   426    103        0 ( 791/ 728- 905/ 770)  18.4 E06-W00 S28-S24
     11 1189  737 (S18W18 L=255)  1115    403        1 ( 984/ 697-1371/ 763)  18.4 W05-W29 S22-S14
     12 1598  629 (S11W47 L=284)   489    136        1 (1544/ 597-1650/ 674)  39.8 W43-W49 S13-S07
     13 1692  623 (S08W55 L=292)   528    108        1 (1669/ 575-1717/ 670)  90.0 W52-W57 S11-S05
     14 1333  564 (S23W33 L=270)  1820    604        0 (1059/ 476-1542/ 622)  14.4 W15-W50 S30-S16
     15 1678  483 (S16W62 L=299)   526     61        0 (1646/ 455-1711/ 507)  59.0 W60-W64 S18-S13
     16 1208  418 (S35W33 L=270)   558    142        0 (1127/ 400-1273/ 439)  29.1 W27-W37 S37-S32
     17 1659  386 (S21W68 L=305)   773    161        1 (1640/ 356-1684/ 419) -14.0 W63-W71 S22-S19
     18 1394  278 (S37W54 L=291)   372    118        0 (1365/ 249-1415/ 309) -15.9 W50-W57 S37-S35
  

    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