Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1597 1692 (N28W71 L= 21)  1242    135        2 (1546/1640-1653/1727) -74.1 W64-W75 N24-N31
      2 1516 1535 (N23W48 L=359)   468     80        2 (1497/1488-1543/1575)  39.8 W44-W51 N20-N25
      3 1366 1544 (N28W38 L=348)   269     83        0 (1328/1518-1401/1569) -59.0 W36-W39 N25-N30
      4 1638 1514 (N18W57 L=  8)   999    175        1 (1563/1466-1711/1547) -45.0 W52-W62 N13-N22
      5  857 1456 (N37W00 L=311)   312     58        2 ( 834/1426- 882/1479) -76.0 W00-W02 N34-N38
      6 1351 1359 (N18W29 L=339)   258     85        0 (1308/1346-1394/1370) -26.6 W25-W31 N16-N19
      7 1767 1165 (S04W52 L=  3)   387     79        2 (1744/1121-1792/1198)  31.0 W49-W55 S06-S02
      8  202 1076 (N26E62 L=248)   587     64        2 ( 183/1046- 225/1112)  45.0 E64-E58 N23-N27
      9  273 1035 (N22E53 L=257)  1913    375        2 ( 174/ 940- 422/1109)   6.1 E66-E38 N17-N26
     10 1577  993 (S08W33 L=343)   327     66        0 (1556/ 968-1595/1028)  36.9 W30-W34 S09-S06
     11 1706  875 (S19W41 L=351)   823    232        1 (1593/ 841-1777/ 915) -17.4 W31-W47 S21-S15
     12 1888  897 (S25W63 L= 13)   619     48        1 (1853/ 878-1920/ 913) -14.0 W57-W69 S26-S22
     13 1827  827 (S27W53 L=  4)   786    109        2 (1793/ 794-1855/ 872)  20.6 W49-W57 S28-S24
     14 1099  838 (S05E00 L=310)   260     50        2 (1075/ 813-1122/ 862) -90.0 E01-W01 S06-S03
     15  522  733 (N01E38 L=272)   528    174        0 ( 501/ 656- 554/ 818)  48.8 E40-E33 S03-N05
     16 1156  550 (S24E04 L=307)   804    219        1 (1094/ 453-1217/ 644)  45.0 E09-W02 S29-S19
     17 1410  463 (S38W12 L=322)   366     78        1 (1389/ 420-1424/ 501)  71.6 W10-W13 S41-S35
     18 1019  263 (S42E24 L=286)  2284    276        5 ( 943/ 176-1074/ 380)  33.1 E37-E14 S49-S33
  

    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