Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1419 1513 (N34W31 L= 89)   435     60        0 (1394/1476-1441/1554) -80.5 W29-W32 N31-N37
      2 1645 1489 (N32W50 L=108)   846    123        1 (1604/1435-1696/1537) -63.4 W46-W52 N27-N35
      3  306 1479 (N32E58 L=360)   922    224        2 ( 273/1444- 359/1517) -15.3 E63-E52 N29-N34
      4  711 1187 (N14E20 L= 38)   248     73        0 ( 680/1157- 750/1216) -36.9 E22-E17 N12-N15
      5 1142 1128 (N10W07 L= 66)   200     74        1 (1106/1116-1185/1144)  11.3 W05-W10 N09-N11
      6 1105  966 (N00W05 L= 63)   283     83        0 (1055/ 945-1153/ 988)   9.5 W01-W07 S01-N01
      7 1849  881 (S08W62 L=121)  1443    285        1 (1815/ 792-1909/ 974)  47.7 W58-W71 S13-S02
      8  526  753 (S13E33 L= 25)   424     96        0 ( 473/ 724- 568/ 781) -23.2 E37-E30 S14-S11
      9  354  622 (S22E50 L=  8)  1468    454        1 ( 272/ 534- 422/ 710)  33.0 E62-E42 S29-S16
     10 1019  575 (S25E01 L= 57)   306     61        1 ( 983/ 559-1052/ 594) -14.0 E03-W01 S26-S23
  

    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