Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  904 1545 (N28E11 L=224)   258    105        1 ( 882/1489- 931/1586) -66.8 E12-E08 N23-N30
      2  776 1486 (N23E19 L=217)   607    167        1 ( 689/1453- 841/1505) -10.3 E24-E13 N21-N24
      3  685 1412 (N18E24 L=212)   225     53        0 ( 651/1389- 713/1430) -21.8 E26-E21 N16-N18
      4  574 1357 (N14E30 L=205)   282     48        0 ( 563/1319- 590/1390) -78.7 E31-E29 N11-N16
      5 1890 1220 (N14W64 L=299)  2064    315        3 (1831/1132-1964/1314) -21.3 W57-W76 N09-N19
      6 1183 1262 (N11W08 L=244)   383    114        1 (1140/1224-1240/1309) -39.8 W05-W12 N08-N13
      7 1564 1216 (N10W33 L=269)   485    155        0 (1497/1196-1622/1250)   0.0 W28-W37 N09-N12
      8 1754 1123 (N06W49 L=284)   321     52        1 (1733/1098-1776/1145) -45.0 W46-W50 N04-N08
      9  432 1040 (S05E38 L=197)   763    151        1 ( 410/ 959- 471/1119)  90.0 E39-E35 S10-N00
     10  829 1056 (S03E12 L=223)   228     52        1 ( 795/1042- 860/1075)   0.0 E13-E09 S04-S02
     11 1852  982 (S01W59 L=295)   396     59        1 (1832/ 950-1872/1008)  45.0 W57-W61 S02-N00
     12  140  905 (S13E67 L=168)  1797    280        7 ( 112/ 815- 170/ 969)  71.6 E71-E62 S18-S09
     13  813  890 (S13E12 L=223)   206     52        1 ( 781/ 874- 845/ 912) -14.0 E14-E10 S14-S12
     14 1113  756 (S20W07 L=242)  1065    251        0 (1015/ 654-1239/ 837)  42.7 W01-W14 S26-S14
     15 1792  552 (S27W63 L=299)  1850    185        1 (1729/ 500-1837/ 625)  48.8 W58-W67 S30-S21
     16  933  547 (S35E04 L=231)   245     47        1 ( 911/ 522- 959/ 570) -14.0 E05-E01 S36-S32
     17  858  498 (S38E09 L=226)   266     89        0 ( 816/ 475- 900/ 530)  33.7 E12-E06 S40-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